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Abstract 


Variations around the average density and composition of the Earth mantle 
may affect long-baseline (anti)neutrino oscillations through matter effects. 
For baselines not exceeding a few thousand km, such effects are known to be 
very small, and can be practically regarded as fractional contributions to the 
theoretical uncertainties. We perturbatively derive compact expressions to 
evaluate such contributions in phenomenologically interesting scenarios with 
three or four neutrinos and a dominant mass scale. 
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I. INTRODUCTION 


The increasing evidence in favor of neutrino flavor oscillations motivates further studies of 
the oscillation phenomenon by means of accelerator neutrino experiments with a (very) long 
baseline L. Current long baseline projects include the Japanese KEK-to-Kamioka (K2K) 
experiment with L = 250 km JjJ, the European CERN Neutrinos to Gran Sasso (CNGS) 
experiment with L = 732 km M (in construction), and the American Main Injector 
Neutrino Oscillation Search (MINOS) experiment with L = 735 km [5J (in construction). 
Possible future projects include oscillation searches with much longer baselines, e.g., using v 
factories |||7J . Although there is still a lively debate on baseline optimization for such future 
projects, there seems to be an increasing consensus in favor of L ~ 2-3 x 10 3 km [|,||. 

For baselines L < 10 4 km, neutrinos cross only the Earth mantle, whose matter density 
can significantly affect the dynamics of flavor oscillations J9fj. For baselines L < 6 x 10 3 
km (as considered in the present work) it has been shown extensively that the main matter 
effects are accurately described by approximating the mantle density profile with its average 
value along the v trajectory (see [|l0|, |n] and the extensive bibliography in |12|). Therefore, 


variations in the matter density (and chemical composition) around the average value can 
be regarded as a fractional contribution to the (theoretical) error budget, rather than as 
a contribution to observable oscillation signals. It should be also stressed that the mantle 
density variations around the average value may also be affected by their own uncertainties. 
Indeed, the standard (radially symmetric) density profile provided by Preliminary Reference 
Earth Model (PREM) jT3| is subject to local variations and uncertainties at the level of ~ 5- 
10%, which may well be nonsymmetrical along a given trajectory [|b| . 


From the above viewpoint, one would like to have a practical recipe to evaluate the 
(small) density variation effects for a variety of possible density profiles (suggested, or at least 
permitted, by local geophysical data), without solving the v evolution equations numerically. 
In this work, by using the perturbative approach applied in (15j for 2z/ oscillations at first 
order in the density variations, we derive such a recipe for 3v and oscillation scenarios 
characterized by a dominant mass scale. We think that our simple expressions may be useful 
to experimentalists involved in simulations of long baseline oscillation signals, in order to 
easily quantify the effects of uncertainties associated to (known or hypothetical) matter 
density variations. 


II. NOTATION AND METHOD 


In this section we introduce the notation and the calculation method. The kinematics is 
defined by the unitary matrix U connecting flavor states u a with massive states v t , 

U* ^ ) U a iVi , (1) 

i 

and by the squared mass matrix 

M 2 = diag {m{, mj ,...) . (2) 

The dynamics is the defined by the hamiltonian 
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M 2 , 

n{x) = u +v{x) 


( 3 ) 


where V is the matrix of potentials felt by neutrinos in matter at the position x (in the 
flavor basis), up to irrelevant terms proportional to the unit matrix 1. 

Solving the dynamics is equivalent to find the evolution operator T from the v production 
point Xi to the detection point Xf, 


Vp{?f) = x i) Va{ x i) , (4) 

a 

which provides the desired flavor oscillation probabilities, 

P(Va -*• Vp) = Pad = \Tp a \ 2 • (5) 


A perturbative evaluation of T can be made 


T5| by splitting the v potential matrix as 


V(x) = V + <5V(a:) , 


( 6 ) 


where 

V— [ f dxV(x) /(xf — Xi) (7) 

JXi ' 

is the average potential matrix along the neutrino trajectory, and SV is the residual variation, 
characterized by 

f 1 dx 5V(x) = 0 . (8) 

JXi 

For later purposes, it is useful to distinguish variations of the neutrino potential which are 
symmetric (+) or antisymmetric (—) with respect to the trajectory midpoint 


x = (xi + x f )/2 , 


namely, 


where 


SV(x) = <5V + 0r) +5V~(a;) , 


5V ± {x — x) — ±5V ± (x — x) . 


The hamiltonian is split as 


H{x) = n + 5V(x) , 


where the constant part is 


— M 2 , - 

H = U ~2E U ' + V - 


( 9 ) 


( 10 ) 


( 11 ) 


( 12 ) 


( 13 ) 
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Correspondingly, the evolution operator is split as 

T = 7 + 5T , (14) 

where T is trivially obtained by exponentiating the constant part of 74, 

T(x f , Xi) = , (15) 

while the correction 8T (at hrst order in hV) is perturbatively given by 

5T(xf, — —if dx T(xf , x) <5V(x) T(x, Xi) + 0(5V 2 ) . (16) 

JXi 

Finally, the oscillation probability is given by 

P<x/3 = |Tpot + ST/3 a | 2 
— Pa/3 + tiPa/3 j 

where 

-P a /3 = \T/3a\ 2 (19) 

refers to the approximation of constant neutrino potential (i.e., constant fermion density), 
while 


(17) 

(18) 


= 2 Re (T.,„ Sl'i 0 ) + 0(SV 2 


( 20 ) 


represents the (hrst-order) correction due to fermion density variations along the neutrino 
trajectory. 

In the following sections we evaluate SP for 3z/ and Av scenarios characterized (from the 
point of view of long baseline experiments) by a dominant (squared) mass scale m 2 , coinci¬ 
dent with the squared mass difference indicated by current atmospheric neutrino experiments 
I6|.[C1| . Our results are valid (at hrst perturbative order) also beyond the approximation of 


a dominant mass scale, provided that the subdominant oscillating terms themselves can be 
treated as a perturbation, which is a good approximation in many cases of phenomenological 
interest for present and future long baseline projects (see, e.g., |18| , [T9|1 ). In this case, the 
effect of matter density variations on subdominant oscillations can be regarded as a second 
order effect, which has been shown to be safely negligible unless L > 7 x 10 3 km |2D], 
case not considered in this work). 


III. RESULTS FOR THREE-NEUTRINO OSCILLATIONS 

Let us consider three neutrino flavors, 

U* = ("e, Vp, v T ) , (21) 

and a mass spectrum with two practically degenerate states, 
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ml ~ ml 


( 22 ) 


plus a third “lone” state which can be either heavier or lighter than the previous “doublet” 
of states, 


2 2 i 2 

m 3 — m 1 2 = ±m . 


(23) 


In the following, the double sign ± is explicitly kept to distinguish the so-called normal 
(+m 2 , upper sign) and inverted (—m 2 , lower sign) spectrum hierarchy. The mass spectrum 
can then be written as 


M 2 


diag 



(24) 


The scale m 2 and the neutrino energy E define the neutrino wavenumber relevant for the 
dominant oscillations, 

k = m 2 /2E . (25) 


The mixing matrix U is parametrized in standard form |22|] as 

u = UM)UM)U 12 {u) , (26) 

where the CP phase S is omitted, being irrelevant under the approximation in Eq. (|2^). 
The potential matrix can be written as 

( V V V\ 

V = diag (+—,——,——) , (27) 

where 


V{x) = y/2 G F N e {x) (28) 

and N e is the electron number density.^ The above notation refers to neutrinos. The case 
of antineutrinos (V —> —V) will be discussed at the end of this section. 

The constant hamiltonian Ti is explicitly diagonalized as 

n = UM)UM) diag (E u E 2 ,E 3 ) UTS)U%iM , (29) 


the energy eigenvalues being given by 


E i 

P/2 

e 3 



=F k — V 



(30) 

(31) 

(32) 


1 In the Earth mantle, it is N e ~ 2 mol/cm 3 . See also Fig. 1 in [|15|]. 
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where 


_ ^20 

■Jicw^V /fc) 2 + s^ 


( 33 ) 


gives the neutrino mixing angle (j) in matter (c = cos, s = sin), and 


k = k s ~it 
S 20 


is the neutrino wavenumber in matter (for average electron density). 
By using Eqs. ([L5|) and ([19|) one gets 


1 Pee Pen Per ^ 


/0 0 0 \ 

Pne P W P n T 

= 1 + 

0-11 

V Pre Prn P tt ) 


o 

1 



(- 1 4 


H- 

4 -4 



\ 4 ~44 


(4 5 23 + 


’20^0 


,2 C2 

5 0 °12 


C 0 

_ q2 2 
*0 C 0 
-r 4 


q 2 _ 02 

*20 °31 


(34) 


(35) 


for the oscillation probabilities with constant (average) potential V, where the oscillation 
factors are given by 


02 • 2 
Sji = sm 



(Xf - Xi) 


(36) 


The above results for P a p are well-known (see, e.g., f23H ). 

We then compute the desired corrections 5P a p due to electron density variations through 
Eq. (|0|). Omitting the algebra, the final results read 


j SP ee 5P eil SP eT \ 
SP,e SP IHl 5P, r = 
V SPre SPt „ 5 P tt ) 


/ 0 0 0 \ 
o-ii 444 cs' 

Vo 1 -l) 

l- 1 4 4 \ 

4 -4 -44 \4m cs " » 
v 4 ~44 ~4 / 


where S' and S" are oscillating terms defined as 


s' = sin y 2 v far - x t )J , 

S" = sin (±|0V ~ x >)'] ' 


(37) 


(38) 

(39) 


while C is basically the Fourier (cosine) transform of the symmetric part of the potential 
variations, 
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C(k) = J f dx 8V + (x) cos (k(x — x)') . 


( 40 ) 


By comparing Eqs. ( |55| ) and it emerges that the effects of nonconstant electron 

density can be simply embedded as corrections to the oscillating factors Sj h 


q 2 
D jl 


Sji + C x 


s|c| S' if jl = 23 or 12 , 
20 


c 9 ^ S" if jl = 13 


(41) 


The above replacement represent our simple recipe for the evaluation of matter density 
variations in the scenario. Once the usual “constant density” probabilities [Eq. (fJ3|)] 
have been computed for a given long baseline experiment configuration, one has just to 
calculate the Fourier term C [Eq. (j4^)] and make the substitutions in Eq. (}4l|) to get the 
corrected probabilities. By using a variety of possible density profiles (allowed or suggested 
by geophysical information along the given baseline) one can then easily evaluate the effect 
of uncertainties in N e upon observable quantities, without solving numerically the neutrino 
evolution equations for the specific profile. 

Several remarks are in order, about the size and the properties of the correction terms 
5P a p in Eq. (|37]). Their size is typically rather small (< 1CT 3 ), both because they are sup¬ 
pressed by an overall factor sin 2 20 (constrained to be ^ 0.1 by reactor^ and atmospheric 
neutrino data [|16[]), and because the integral C is typically rather small, due to the oscillating 
behavior of the integrand. The maximum of C is reached only for a hypothetical config¬ 
uration of the matter prohle with 5N+(x) oc cos[/c(a: — 3c)], namely, when matter density 
fluctuations happen to be “in phase” with oscillations. In such a contrived case, since 


dx8V + = 0.386 x 10 



(42) 


with (plausibly) 5N e < 10%N e — 0.2 mol/cm 3 , the integral C could be as large as a few 
percent over long baselines (> 10 3 km). However, the realization of such a “maximum 
effect” configuration (if approximately reachable) can be attained only in a narrow range 
of k (i.e., of neutrino energy). Outside such a range, the size of C (and thus of 5P a p) is 
rapidly suppressed. We refer the reader to the recent papers |^4[[T2j for extensive numerical 
investigations of the size of 5P a p, for various admissible SN e (x ) profiles. 

Concerning the symmetry properties of the oscillation probabilities, notice that the first- 
order corrections 8P a p do not involve the antisymmetric part of the potential. As a conse¬ 
quence, T-violation effects are absent not only at 0th order, but also at 1st order in 5V, 


P a(3 P Pa j 

Pa 0 = Ppa + 0(SV 2 ) , 


in agreement with the discussion in [2C ,21 


(43) 

(44) 


2 Although reactor bounds apply, strictly speaking, to <f> and not to q i>, our arguments are not 
qualitatively changed. 
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Finally, let us consider the case of antineutrinos , for which the potential is given by 
—V. It is easy to check that all the above oscillation probabilities are invariant under the 
simultaneous replacements ±m 2 —> =pn 2 (change of hierarchy) and +V —> —V (change of 
potential). Therefore, antineutrino probabilities for direct (inverse) hierarchy are equal to 
neutrino probabilities for inverse (direct) hierarchy, 

P(±m 2 ; V a -> Vp) = P(=Fm 2 ; vp -> u a ) . (45) 


IV. RESULTS FOR FOUR-NEUTRINO OSCILLATIONS 


In this section we consider the phenomenologically interesting scenario of three active 
plus one sterile neutrino u s , with a mass spectrum consisting of two separated doublets (see 
and references therein). In such a case, it is useful to decompose the flavor space as 

(z/ e , Up, u T , v s ) = (i/+, © (z/_, v e ) (46) 

where the states v± are linear combinations of u T and u s , 



V- 


c e V \ 

-S ? C* ) \U S ) 


(47) 


Under plausible approximations discussed in [17[], the mass-mixing parameter space can then 
be factorized as 


4/z parameter space ~ (m 2 , V 1 , 0 ® (5m 2 , u>, £) , 


(48) 


where the hrst and second factors dominate the oscillation physics in terrestrial and solar 
experiments, respectively. From the point of view of long-baseline experiments, the role of 
(m 2 ,^) is similar to the 3u case in the previous section, while the new angle £ interpolates 
smoothly from pure active oscillations (s^ = 0) to pure sterile oscillations (sg = 1). 

In the (v + , Vp) basis, the effective hamiltonian for long baseline experiments is given by 

n 

^ = 1 ( v - k C 2 y _k S 2 y 
2 \ ks 2 ^ -V + kc -2 v, 



where the average neutrino potential is now related to the average neutron density N n , 

V= V2G F s\N n /2 . (50) 

The dynamics is then effectively reduced to an equivalent two-flavor problem, which is 
easily solved. We give the final results for the relevant muon neutrino probabilities, 


-P/4/4 — 1 


s 2 - S' 2 

b 2ij) ° 


P — c 2 s 2 - S' 2 

-T/ 4 T — ° > 

~p _ „2 2 _ o2 

r /4S — b 2ip ° ) 


(51) 

(52) 

(53) 












where 


s 2\p 


and the dominant oscillating factor is 


S2 v> 


(c 2 „ - V/kf + S 2 


(54) 


2i/> 


S 2 = sin 2 (-(i/-li) I , 


(55) 


with 


k = k -^ 

s 2 tji 


(56) 


We perturbatively find that the first-order correction due to nonconstant neutron density 
amounts to the following replacement 


s 2 -* s 2 +c 2 fC sin (hi/ -Xi )1 , 


(57) 


where C is formally defined as in Eq. (|40|). 

Notice that, contrary to the previous 3u scenario, the cases +m 2 and — m 2 are now 
equivalent, being connected by the replacement 0 —> 7r/2 — -0. One can then fix the sign of 
m 2 , provided that 0 is taken in its full range [0,7r/2]. 

Finally, antineutrino probabilities are obtained by the above neutrino probabilities 
through the replacement V —► —V or, equivalently, through m 2 —> —m 2 or, equivalently, 
through s^, —> Cl/,. 


V. SUMMARY 


In the context of long baseline v oscillation experiments, the small effects due to non¬ 
constant matter density (along a given neutrino trajectory) can be regarded as a fractional 
contribution to the uncertanties associated to oscillation signals. Therefore, it makes sense to 
search for a practical method to compute them for given density profiles. We perturbatively 
find that the first-order effects can be simply embedded as a correction to the oscillating 


factors in both 3 v and Av scenarios with a dominant mass scale [Eqs. (fflD and (|57|), re¬ 
spectively], The correction involves the Fourier (cosine) transform of the symmetric part 
of the density variations, which can be easily evaluated and incorporated in experimental 
simulations. 
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